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Abstract 

We address the issue why a cosmological constant (dark energy) pos- 
sesses a small positive value instead of being zero. Motivated by the 
cosmic landscape picture, we mimic the dark energy by a scalar field 
with potential wells and show that other degrees of freedom interact- 
ing with it can localize this field by decoherence in one of the wells. 
Dark energy can then acquire a small positive value. We also show 
that the additional degrees of freedom enhance the tunneling rate be- 
tween the wells. The consideration is performed in detail for the case 
of two wells and then extended to a large number of wells. 



1 Introduction 



Observations indicate that our Universe is currently accelerating [T]. The 
simplest explanation of this fact in terms of the Einstein equations of gravity 
is the existence of a non-zero Einstein cosmological constant A. However, a 
more general dynamical entity dubbed "dark energy" (DE), whose effective 
energy-momentum tensor is close to that of A, is also a serious possibility 
[21 El IH El [7] . Present measurements of the effective DE equation of state 
w~de '■= Pde/pde are inconclusive about its nature and are compatible with 
DE being given exactly by A, that is, with w DE = — 1 [TJ. 

The microscopic origin of DE is presently unknown and can perhaps only 
be explained after a full quantum theory of gravity is available [S]. There 
are actually two aspects connected with DE. First, it is unclear why the 
effective energy density of DE (or A itself) is much smaller than one would 
expect from elementary particle physics on dimensional grounds. Second, it 
is unclear why it is not exactly zero but has a small positive value. The first 
question might only be answered on the basis of a full quantum theory of all 
quantum fields including gravity. But it is possible that the second question 
can be addressed on the basis of known physics. 

A first step in the understanding of the second question was made by 
Yokoyama [5J. He modelled DE by a scalar field <fi with a potential V{<p) 
that is characterized by a double well. This choice can be motivated by re- 
cent ideas in string theory, where a "landscape" of many (perhaps as many 
as 10 500 or more) local minima of a complicated potential is discussed, see, 
for example, [101 [11] and the references therein. Yokoyama assumed that, 
perhaps due to some unknown symmetry, the exact ground state of the Uni- 
verse is characterized by a vanishing vacuum energy, that is, a vanishing 
cosmological constant and that the observed small deviation from zero could 
arise from the fact that the Universe is not in its ground state. 

But how can this happen? The ground state for a double-well potential is 
extended (delocalized) over both minima. In contrast to this, a state localized 
in one of the minima is a superposition of the eigenstates; in the simplest case, 
it is a superposition of the ground and the first excited state. The effective 
energy of such localized states is greater than the ground-state energy and 
would therefore be positive if the ground-state energy were zero. If the wall 
between the wells is not too small, the values for this positive-energy states 
are tiny because they differ from the ground-state energy only by a small 
tunneling factor proportional to exp(— So), where Sq is the instanton action. 
The reason for the observed small positive cosmological constant could thus 
lie in the fact that the dark energy in the Universe is in a localized state that 
is concentrated near one of the minima of the potential. 
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As long as the universe stays in such a localized state, the effective DE 
equation of state will be u>de ~ — 1- There exists, however, a certain proba- 
bility that the Universe can evolve into the ground state, which is a super- 
position of the localized states, or tunnel into another localized state. The 
question then arises how large the time scale and the tunneling rate would 
be. Obviously, these scales should obey all known observational constraints 
presented in [I]. 

In our paper we shall elaborate on this idea in two respects. First, it 
has to be justified why a DE field is not found in its ground state in the 
first place, but in a localized state. This is reminiscent of an analogous sit- 
uation in quantum mechanics where chiral molecules are typically found in 
left-handed or right-handed isomeric forms and not in a superposition of the 
two (as e.g. the ground state would be). Historically, this was called "Hund's 
paradox" [12]. Its solution is based on the central concept of decoherence 
[T2l [13] and was recently presented for a specific case in quantitative detail in 
|14j . Decoherence is the unavoidable and irreversible emergence of classical 
properties with its environment; "environment" stands here for any irrele- 
vant degrees of freedom that interact with the quantum system and thereby 
become entangled with it. In the case of the chiral molecules these can be 
photons or air molecules scattering off these molecules. Decoherence will also 
play a central role in our analysis, in order to explain why the DE field is not 
in its delocalized ground state, but in a localized state with higher energy. 

Our second elaboration is a direct consequence of the first one. If addi- 
tional "decohering" degrees of freedom are present, they will have an effect 
on the tunneling rate from one well to the other. We shall thus discuss both 
the pure tunneling rate of the isolated system as well as its modification 
by the environment. In view of the experience from quantum mechanical 
models [T5], one would expect that these degrees of freedom will in general 
reduce this rate, so that tunneling will become less likely. In contrast to this 
expectation, however, we shall find that this tunneling rate is here increased 
rather than decreased, by a field-theoretic mechanism similar to the Casimir 
effect. The localization by the environment together with the magnitude of 
the tunneling rate could then provide the explanation of why we observe a 
small positive value for dark energy. 

We note that our approach is completely different from the one used in 
the paper [16] (which appeared when our paper was prepared for publica- 
tion), where entanglement between different cosmological epochs, not differ- 
ent vacua, is proposed as a source of A. 

Our considerations should also be relevant for the inflationary stage in the 
early Universe, which was dominated by a primordial DE whose properties 
were close to an effective (though metastable) cosmological constant, too. 



2 



However, we do not attempt here an application in this direction. 

Our paper is organized as follows. In Sec. 2.1 we present our model: 
DE is modelled by a scalar field (quintessence), and it interacts with an 
environmental massless scalar field that may arise from metric perturbations. 
In Sec. 2.2 the process of decoherence is discussed in detail for both sub- 
and super-Hubble modes of the environment. It is found that these modes 
localize the DE field in one well, the super-Hubble modes even more efficiently 
than the sub-Hubble modes. In Sec. 3 we calculate the tunneling rate for a 
Minkowski background as well as an expanding Universe. We find that the 
presence of an environmental field enhances the original tunneling rate. In 
order to approach more closely the ideas of a cosmic landscape, we extend 
our discussion in Sec. 4 to the presence of many vacua. A brief Appendix 
summarizes some material about the (^-function renormalization needed in 
Sec. 3. 

Our main results have already been briefly announced in [T7j . 

2 Small positive dark energy from decoher- 
ence 

2.1 The model 

We present here our model in which a scalar field mimicking DE (called 
the "system" in the following) is coupled to other degrees of freedom called 
"environment". As in the Yokoyama paper [9J, 4> is assumed to possess 
a potential with two quasi-localized minima. As for the environment, we 
choose another scalar field called a which couples to <fi. The details of the 
interaction are of secondary importance; it only has to be able to generate a 
sufficient entanglement between system and environment [12]. 

Both system and environment are supposed to evolve in a flat Friedmann- 
Robert son- Walker (FRW) background with the line element 



{h — c — 1 is set throughout the paper). We assume here for later convenience 
that the scale factor, a, has the dimension of a length, while x, y, and z are 
dimensionless. The total action of system and environment then reads 



where g is the determinant of the metric g^ v , and £ sys , C cnv , and C mt denote 
the Lagrangian of the system, environment, and interaction, respectively. 



ds 2 = g^dx^dx" = dt 2 - a 2 (t)(dx 2 + dy 2 + dz 2 ) 



(1) 
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The spatially homogeneous scalar field describing the vacuum energy reads 

£ sys := \<P 2 ~ V(<j>) , (3) 

while the Lagrangian of the environment, £ env , describes a massless scalar 
field cr, 

C m := l -g B ^d^d v a. (4) 

The interaction between system and environment must be able to dis- 
criminate between different values of the scalar field (p. In quantum mechan- 
ical applications one often chooses a bilinear interaction, see for instance the 
Caldeira-Leggett model [J5] or the spin-boson model described, for example, 
in Section 5.3 of [13]. Since, however, our field <fi is spatially homogeneous, 
a bilinear term of the form J d 3 x er(x, t)<p(t) would give no interesting dy- 
namics, since only a single Fourier component of the scalar field a, the one 
with vanishing momenta, would interact with the system field. We therefore 
introduce the tri-linear interaction 

Cmt,a ■= -g s o- 2 (p (5) 

for the coupling with the scalar-field environment. Here, the coupling con- 
stant g s has to be chosen such that the product g s cj) is positive. This guar- 
antees that the corresponding Hamiltonian is bounded from below and that 
the dynamics is thus stable. Since both a 2 and </> 2 have physical dimension 
mass (M) over length (L), the dimension of g s is (ML 3 ) -1 / 2 , which in the 
natural units used here is equal to M. 

What could be the origin of such a coupling? It can arise, for example, 
from the expansion of the metric determinant \f—g into the scalar and tensor 
modes. Let us consider therefore the FRW line element with scalar and tensor 
perturbations (see e.g. [IB]). 

ds 2 = (1 + 2ip 1 )dt 2 - a 2 [(1 - 2ip2)Sij + Kj] dx { dx j , (6) 

where ^1/2 are scalar perturbations and fay are tensor modes. In the trans- 
verse and traceless gauge, there are only two independent tensor modes, 
^23 = ^32 and /122 = — ^33- If there is a linear term in the potential of the 
system field 0, we get from the expansion of the determinant a term of the 
form 

J d 3 xy/^g<p « (7) 
J d 3 xy^ (l + ^ - 3^2 - 2^1 + - ^2? - \h 2 22 - • 
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Discarding the terms linear in ipx and ip2 for the reason mentioned above, we 
are left with tri-linear interactions of the form §5§ ; the role of our field a could 
thus be played by ifti, ip2, ^22, or /i 23 . Such interactions can thus arise both 
from the scalar and the tensor perturbations of the metric itself, although it 
is conceivable that a is any other field that occurs in a fundamental theory. 

We shall now address the Hamiltonian that will be used in the quantum 
theory. It can be derived from (121)— (15]) and in the momentum representation 
reads 

-P) + ^(p)a(-p) (8) 

+ , 
where 

*(x) = J rfW ab (pK xp , (9) 

n(x) = ^|^n(p)e^, (io) 

and denotes the pure 0-Hamiltonian (see (fTTj) below). Note that p is 
dimensionless because x is dimensionless. 

In the following, we shall simplify the part of the Hamiltonian describing 
the system. We shall assume that the dynamics of the system is dominated 
by the two lowest energy eigenvalues of the double-well potential, that is, we 
assume that their difference is much smaller than the energy gaps within a 
single well. It is then possible to reduce the system to an effective two-state 
system, 

where E + and E_ are the energy levels of the localized minima, and A is the 
tunneling matrix element. 

The reduction of the system to an effective two-state system leads to the 
interaction 

iW = g.a*{t) ( O + M J d 3 p a(pM-p) , (12) 

where the environment can only discriminate between the two different min- 
ima of the potential, <p + and </>_. 
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= I dp A \ ^n( P )n 



+# s a 3 0a( P )a(-p) 
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All together, our model resembles a spin-boson model [13], although the 
coupling in the standard situation is taken to be linear in the environmental 
fields. It is well known that situations with a double-well potential can often 
be described by an effective two-state system [TJl 03] • I n this context the 
Hamiltonian (llip would be written in the form 

= a x A + -(SE)a z + -(E + + E.)I , 

where a x and a z are Pauli matrices, and SE := E + — E_. Spin-boson models 
describe the interaction of a central system with its environment in the case 
when the system is effectively acting as a two-level system. 



2.2 The reduced density matrix 

With this simplification at hand, it is possible to calculate the reduced density 
matrix of the two-state system. We assume that the initial state is a product 
of a system and an environmental state, |^) sy s <8> |^) C nv The time evolution 
will then generate an entanglement between them. This evolution is governed 
by the functional Schrodinger equation 

= H^\V) . (13) 



In this section, we shall neglect tunneling, that is, we set A = in (JUj) ; this 
enables us to solve the Schrodinger equation exactly. We assume that the 
state of system and environment is of Gaussian form (see e.g. [191 1201 12T] ) 
and make the ansatz 

faN + (t) exp {-\jd 3 p a(p)n+((f>,p, t)a(-p) - iE+t)\ 

m=[ , (14) 

\(3N_(t)exp {-Ud 3 p a(p)n_((f),p,t)a(-p) - iE.t) ) 

where Q + /_(p, t) and N + /_(t) are time-dependent functions to be determined 
from the Schrodinger equation, and a and /3 are constants with |a| 2 +|/3| 2 = 1. 
With the above ansatz one obtains the following Riccati-type equations, cf. 



(fi +/ _(p,t)) 2 

73 



-/-KP^JJ 2 3 

- L) = 

and 



(p, t) = ^ 3 + p 2 a + g s a 6 <f) +/ - (15) 



^TVSl^tp,^^. (16) 
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Eq. ( 1T5|) can be transformed by the ansatz 

fi +/ _(p,t)=:-m 3 ^M (17) 

into a linear equation for «+/_. Switching to conformal time, which is defined 
by a dr] = dt, and denoting derivatives with respect to t] by a prime, we obtain 

2a' 

2 + 7a 2 )u +/ _0, 77) + —«'+/_ (p, 77) + u"+/_(p, ^) = , (18) 
where we have introduced the quantity 

7 := 2^ s +/ _ , (19) 

which has dimension L~ 2 = M 2 and obeys 7 > 0. (We have skipped the 
indices +/— in 7 for simplicity.) We assume that the coupling between 
system and environment is small, that is, 7 < H 2 . The total density matrix 
for system and environment is given by the pure state 

p(f) = |*)<*|, (20) 

from which the reduced density matrix is obtained by integrating out the 
environmental scalar field, 

p{t) sys :=Tr onv |tf>(#|. (21) 

In position representation, this reads 

Pij = Jva , (22) 

where i,j run over the values + and — , and \& = ip-) T - Since by setting 
A = we have neglected dissipation, the diagonal elements of the reduced 
density matrix remain unchanged, that is, one has 

p++ (t) = |«| 2 |iV + (t)| 2 J Va exp - J d 3 p a(p)Ml+(t)a(-p) 

= \a\ 2 = P++ (0) (23) 

and analogously p (t) = p (0). The probabilities of finding the system in 

state + or — are thus unchanged by the environment; this corresponds to a 
quantum-nondemolition (or ideal) measurement 
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The non-diagonal elements can be calculated as follows: 
P+-{t) = P *_Jt)=aP*N + (N_)*x 



Daexp 



d 3 p tr(p) (n+ + cr(-p) - *(£+ - E_)t 



a (3 



det 1/4 (Kn+)det 1/4 (3?fi_ 



x 



det 1/2 ((ft+ + (ft„)*)/2) 
xexp^-^y ^ / ^Tr(SRn+ - 3fft}_) - z(£+ - #-)*V (24) 

From (fT8|) we can see that the functions f2 + /_ depend on the small parameter 
a 2r y/p 2 for sub-Hubble modes and on j/H 2 for super-Hubble modes, see 
below for an explanation of these terms. Expanding f2+/_ up to second 
order in 7 yields 



n+/-(7) 



+/- 



7=0 



d o 



7=0 



2^ 



7 



7=0 



-: f2 + f2 7 7 + -f2 77 7 2 



(25) 



where derivatives with respect to 7 are denoted by indices. (Strictly speaking, 
the expansion is with respect to the dimensionless combinations a 2r y/p 2 or 
l/H 2 .) 

The approximate expression for the non-diagonal elements then reads 



P+-{t) = P+-(0)exp 
with p + _(0) = a/3* and 



-Tr 



(m,) 2 



Tr 



v Kfi 2 

' \ 1 
(it — Tr(KQ_ 



+ 



30 



77 



m.) + (^ 



(my 

(Kfi) 2 
-#_)t. 



(26) 



(27) 



In the following, we want to discuss the explicit form of the decoherence 
factor. To begin with, we consider the impact of the sub-Hubble modes on 
the system, that is, the impact of modes whose wavelength A = 2ira/p is 
smaller than the Hubble scale H~ x . Using a WKB-approximation, which is 
adequate for p 2 3> a" /a, the solutions to the differential equation (ITS]) read 



u 



+/- 



A ( P 
- (p, v) = —expli dfj u +h (77) 



(28) 



where 



//\ 1/2 




"+/-(77):= (P' + T« --) • (29) 

We have chosen an adiabatic vacuum in the infinite past for each Fourier 
mode p. This choice is possible because the interaction vanishes for a — > 
when the modes are far inside the horizon. In the case of a massive scalar 
field (g s (j) + /_ > 0) in the exact de Sitter background, it leads to the Bunch- 
Davies vacuum for the whole field a, see e.g. [22] • The trace of the real part 
of the exponent in (1261) is 

(30) 



where V denotes the dimensionless coordinate volume which must be fixed 
by an appropriate infrared cutoff. The WKB-approximation holds for modes 
far inside the horizon, p min > yfa"/a\ in the case of a constant Hubble rate, 
we have p m i n > \/2Ha. 

To discuss the explicit form of the decoherence rate for super-Hubble 
modes, that is, for modes with wavelengths greater than the Hubble scale, we 
shall restrict ourselves to the expanding de Sitter case where H = constant. 
In this case, it is not possible to find WKB solutions valid for all times in- 
cluding t — > oo (77 — >■ 0) , because the time evolution for super-Hubble modes 
is highly nonadiabatic. The solution of Eq. (ITS]) that has the correct WKB 
behavior for rj — > — 00 is then given by 

« +/ _fe,) = ^g) 3/2 W« 7JI ^W, (31) 

where denotes a Hankel function. (Recall that 7 < H 2 .) In the massless 
case g B — 0, this expression reduces to the textbook free field solution for the 
adiabatic (WKB) initial vacuum state at rj — > — 00: 

^^=-^{ i+ £) eipv - < 32) 

Note that due to the absence of WKB solutions for t — > 00 (rj — > 0), there 
exists an ambiguity in the definition of a particle for super-Hubble modes for 
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sufficiently small masses m 2 = 7 2 < 9H 2 /A. In addition, these modes are 
strongly squeezed. In particular, in the massless case g s = we have 

2 a 2 H 2 

smh r p = ' ( 33 ) 



where r p is the squeezing parameter of the mode, see e.g. [231 EI]- On the 
other hand, in the exact de Sitter case it is possible, in principle, to consider 
the whole solution (1321) as the second-order corrected WKB solution that 
would imply no massless "particle creation" in the de Sitter background, see 
e.g. (23]. However, any ambiguity disappears if we take into account that 
H is not constant in any realistic inflationary model and decreases after the 
end of inflation (or even during it), see the more extended discussion of this 
problem in [26J. 

Using ( 133]) . we can express f2(p, 77) entirely as a function of the squeezing 
parameter: 

"<"■'> =i=5S£v (34) 

Using ( 153|) and the approximation of (T5T]) for pr) — > 0, 

" + /-(™)-g) 3/2 (f) _V ^. (35) 

we obtain for the impact of super-Hubble modes on the system the result 
„ /(^ 7 ) 2 + (^ 7 ) 2 \ Va 2 f*-*, ^ , , . , 2 , 2 

Ua 2 



2a 2 H 2 2a 2 H 2 \ ( a 4 H 4 a 4 H 4 



Pmax Pmin V V^Pmax 3p min 



(36) 



This result for the trace depends on the minimal and maximal values for 
the dimensionless wave number. For the super-Hubble modes, we take for 
the minimal wavelength the Hubble scale, so p max = 2iraH. What has to 
be taken for the maximal wavelength? Clearly, an infinite wavelength would 
lead to p min = and thus to a divergence in fl35J) . In a closed universe, a 
reasonable cutoff would be A = a. In the case of an open universe, we shall 
adopt the argument on p. 159 of [27], which goes as follows. We assume that 
the initial fluctuation spectrum has a cutoff at Aq ~ H^ 1 , where the initial 
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size of the inflating region is about H^ 1 . This leads to a cutoff A = a as in 
the closed case and we can set p min = 2n in (156% 

However, in the case of a real post-inflationary universe, it is natural to 
take for the cutoff of A its homogeneity scale, that is, the scale, at which per- 
turbations generated during inflation become of order unity and the space- 
time description using the FRW background (jTJ) loses sense. Though this 
scale is much smaller than the radius of pre-inflationary curvature, it is still 
much (typically exponentially) larger than the Hubble scale after the end of 
inflation. 

Evaluating the trace using these numbers and inserting the result into 
(126|) . we find for the absolute value of the non-diagonal element of the density 
matrix: 



\P+-{t)\ = |p+-(0)|exp 



72n 2 H 2 



24tt 



0(aH)' 



• (37) 



We recognize explicitly that this non-diagonal element becomes increasingly 
small for increasing aH, that is, decoherence becomes efficient and the field 
is localized in one or the other well. 

Evaluating in the same manner the density matrix (1301) for the sub-Hubble 
modes, we obtain instead 

= |P+-(0)|exp I q^2jj J ' ^ ^ 

where C ~ 0.329575. Taking the ratio of the widths of the two Gaussians 
(I37jl and (1551). we get 



(A0) s 2 upcr 8.97T 3 
(A0)L * W ' 



(39) 



which goes to zero for aH — > oo, that is, the super-Hubble modes are much 
more efficient in the localization of than the sub-Hubble modes. The reason 
for this is that the ensuing entanglement with the DE field is stronger with 
the super-Hubble modes due to the stronger interaction. 

One can associate with fl37|) a decoherence time roughly as follows. As- 
suming that a(t) = H^ 1 exp(Ht), the condition that the (absolute value 
of the) exponent becomes of order unity leads to a decoherence time ~ 
(6H)~ 1 times logarithmic terms containing the coupling and the separation 
of the minima (counted from the beginning of the de Sitter stage). 

The environmental field a used in the investigation above can have differ- 
ent origins. The most natural one is given by the scalar and tensorial modes 
of the metric itself. The dynamics of these modes during an exponential 
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expansion of the background is well known and can be described, for exam- 
ple, in the squeezed-state formalism already mentioned above [23J. If the 
dynamics of these modes is unitary, they become highly squeezed in the field 
momentum and highly elongated in the field variable. However, these modes 
are themselves prone to decoherence [28j EH [26] . The classical pointer basis 
distinguished by the interaction of these modes to their environment is the 
field-amplitude basis; this is why one can describe the primordial fluctuations 
by classical stochastic field variables. Why, then, can these fluctuations serve 
as a quantum environment to decohere the DE field? 

Let us assume that the modes a p couple to another field with modes Xk 
responsible for their decoherence. Unless the coupling is very strong (which 
would not be realistic), the total quantum state is of the form 

M, { X k}) ~ M<l>, W P })MM, {xk}) . (40) 

The entanglement in the wave function ^({cp}, {Xk}) is responsible for the 
decoherence of the modes a p [28]. However, tracing out the modes Xk i n the 
full state (I4UI) is ineffective as far as the DE field is concerned; for the deco- 
herence of the latter only the wave function ipi(<f), {c p }) is responsible. That 
is, our results presented in this section are insensitive of the modes a p being 
themselves decohered or not (provided, of course, the involved interactions 
are not too strong). 

This is in accordance with decoherence in quantum mechanics [T2]. Scat- 
tering processes can localize a quantum system such as an electron or a 
molecule. The important point is that the scattering causes entanglement 
between the scattered system and the scattering agency and that the infor- 
mation about the superposition is no longer available at the system itself. 
Typically, the environment is itself decohered, but since the full states are 
usually of the form (!4Up . this decoherence is irrelevant for the decoherence of 
the original system. A somewhat related example is the one discussed in [29] . 
There, two oscillators become entangled with the same heat bath. However, 
unless the two oscillators are close together, the respective entanglements 
with the bath will not lead to an entanglement between the oscillators them- 
selves, that is, the total state will be of a form similar to ( l4"Uj) . 

To summarize, we have shown in this section that it is justified to assume 
that the DE field is localized in one of the two wells only; interference terms 
between the two wells are dynamically suppressed by decoherence. This justi- 
fies the scenario presented by Yokoyama that the small positive cosmological 
constant could arise from the DE field not being in its ground state - it is 
localized by decoherence in one of the minima of the potential. 

An extension of Yokoyama's work to the case of many wells, taking into 
account ideas from string theory, was suggested in [30]. There, the authors 



12 



assume the ground state of our Universe to be a superposition of all ac- 
cessible vacuum states. However, this seems to be a doubtful assumption, 
since the unavoidable interaction with environmental degrees of freedom, for 
example Standard Model fields or thermal excitations, should lead to their 
decoherence. 

3 Tunneling rate 

In this section we shall investigate the influence of the system-environment 
interaction on the tunneling rate. Before we discuss this in detail, we want 
to recall some basic facts about tunneling in field theory. We shall first 
address the case of a Minkowski background and then turn to the expanding 
Universe. 

Following [31], we know that the tunneling rate of a system given by a 
scalar field Lagrangian is of the form 

T = Aexp(-S E ( ( p)) , (41) 

where Se{4>) is the classical Euclidean action of the scalar field evaluated 
along the tunneling trajectory of <p\ the prefactor A can be determined by 
the second variation of the action. 

According to the classical equations of motion, the field adopts for most 
of the time the value cj) + of the false vacuum and approaches the value 0_ of 
the true vacuum after a short transition time. The terms "false vacuum" and 
"true vacuum" may be misleading, since they denote the classical minima of 
the potential (with + representing here the higher minimum), in contrast 
to the true quantum mechanical vacuum which is a superposition of 0+ and 

0- 

The tunneling time T between the two vacuum states is assumed to be 
large compared to the characteristic instanton transition time 1/u. We thus 
consider situations in which various tunneling processes from one minimum to 
the other can be considered separately. This picture of separated transitions 
can only be justified by decoherence, since the superposition principle is 
universally valid and thus holds also for widely separated "jumps". 

Assuming spherical symmetry, a transition between the localized vacuum 
states can be described by the growth of a bubble that is nucleating sponta- 
neously at a radius R . The true vacuum inside and the false vacuum outside 
the vacuum bubble are separated by a wall with a negative surface tension. 
In the limit of a small energy difference between the localized vacuum states, 
the nucleation radius is given in terms of the energy difference and the surface 
tension [31]. 
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The interaction with the environment will influence the decay rate due to 
decoherence and dissipation. The authors of (T5] considered a macroscopic 
position variable coupled to a heat bath of harmonic oscillators. After inte- 
grating out the environmental degrees of freedom, they obtained a Langevin 
equation with an effective friction term resulting from dissipative effects. The 
consequence of this friction term is to reduce the tunneling amplitude. 

In contrast to [T3] we consider the tri-linear interaction (jSJ) between sys- 
tem and environment. In addition we work in the context of field theory and 
thus do not restrict ourselves to a quantum mechanical model, that is, to a 
finite number of oscillators. 

Starting from (jSJ) and switching to Euclidean time, t — > —it E , we find the 
Euclidean action 

S E = S E , sya + J dx 3 dt E a(x, t) y~0 E + g s <f>(t E )j <r(x, t) . (42) 

Integrating out the environmental field a leads to the following formal ex- 
pression that modifies the tunneling rate fj4TT) . 



r = T N x / ZMx, ( )exp(- SE ) = ,/_^zS£L_ ro . (43) 

The normalization TV was chosen such that r| Ss=0 = T . Solving the nonlocal 
equations of motion given by the variation of 

S E ,eS = S E)Sys + -Trm(-D E + 2g s cf>) (44) 

and evaluating the effective action along the tunneling trajectory of <fi would 
give the exact modified tunneling amplitude. In order to simplify the cal- 
culations, we want to assume that at lowest order the trajectory of 4>{t E ) is 
given by the unperturbed equations of motion. In this approximation we can 
compute the functional determinant. 

For this computation we have to solve the eigenvalue equation 

~-^ + 2g 8 <f>{t E ^ = \tl). (45) 

Since the nucleation of the bubble is spherically symmetric, it is appro- 
priate to use the Laplacian in spherical coordinates. Making the ansatz 
= 4 > in(r)yim{0,(( ) )u(t E ), we find for the eigenvalue equation of the radial 
component 

■^!T r 7T + 2 <t>in{r) = K nm <j) ln {r) . (46) 
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A natural boundary condition would be 4>i n (Ro) = 0, that is, the eigenfunc- 
tions are vanishing at the boundary of the bubble. Equation 046|) is solved 
by the spherical Bessel functions, that is, <f>i n (r) oc ji(Ki n r). The eigenvalues 
Kin are the n-th root of ji divided by Rq- We thus have 

°° / d 2 \ 

Det(-D E + 2g s (j){t E )) = J] {21 + l)Det nZ ( + n 2 n + 2g a <j>{t E ) J , (47) 

n,l=0 \ E J 

where the degeneracy of the eigenvalues was taken into account with a factor 
2Z + 1. 

Since the eigenvalues of the spherical Bessel functions are not explicitly 
known, we simplify the problem by assuming periodic boundary conditions 
in a volume L 3 with L = 0(Rq). For the spatial part of the Euclidean 
d'Alembert operator, we choose periodic boundary conditions with a length 
L. The functional determinant involved in (j43j) then separates into a product 
of functional determinants labeled by the mode number n: 

Det(-D E + 2g s <f){t E )) = _Q 47rn 2 Det n f + + 2g s <f){t E ) J . (48) 

n=0 V S / 

This expression is divergent for two reasons. First, for each fixed mode 
number n, the determinant is an infinite product of eigenvalues, which is in 
general infinite. Second, due to infinitely many modes, the situation becomes 
even worse. 

In order to regularize the expression ( 148|) we choose the ^-function reg- 
ularization method presented in [32]. A short summary of this method is 
given in the Appendix. For any second-order differential operator T> we can 
write 

(DetX?) 1 / 2 = exp (-^(0) " ^C(O) hi/i 2 ) , (49) 
where £(s) is the generalized zeta function 

coo = E^ ( 5 °) 

involving all eigenvalues A of the differential operator. The parameter \i 
appearing in ( 149]) is a renormalization parameter with dimension of a mass. 
According to [32], the £- function reads 

t_ , . sin(Trs) f°° dM 2 d 2 T .^ 2 . 
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with 



I(M 2 ,s) = J2 — lnu(M 2 n 2 ,t E ). 



(52) 



n=l 



The integral flBTl) converges for some s > 0, since I(M 2 ,s) increases with 
finite polynomial order proportional to M k , and can be analytically continued 
to s = [32]. 

The functions u(— A, are the eigenf unctions of the differential operator 
under consideration. The eigenvalue equation corresponding to (]4"gj) reads 



Note that the differential operator is always positive definite since g s 4>{tE) > 0. 

In general one needs two independent boundary conditions in order to 
determine the eigenvalues of (153"]) uniquely. We assume that the mode func- 
tions have a root at nucleation time To, that is, u{— A, To) = [32]. This 
time is usually of the order of or equal to the nucleation radius Rq [31] . The 
second boundary condition is given by the normalization, see below. Since 
the regularization method employed discards one of the two independent so- 
lutions of equation (153"]) (see the Appendix for a detailed discussion), the 
eigenf unctions are uniquely determined by a normalization condition. 

The leading term of the uniform WKB-solution of (153]) reads 



where the choice of tE,o determines the normalization of u{— A, We as- 
sume that the scalar field <p is in the true vacuum at some large positive time 
tE, that is, (j)(t E ) = <p-i and fix the normalization such that 



Since the normalization of u{— A,£e) enters the function J(M 2 , s) only log- 
arithmically, the choice of a different normalization will not significantly 
change the results. 




(53) 




(54) 




(55) 
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For large mode numbers n there is an approximate degeneracy of 47m 2 . 
We find 



I{M\s) 



^ 47m 2 



n=l 



n 



2.s 



2^2 



47T n 
_ L 2 " 



+ 2£ s 0(t s ) + M 2 n 2 



47r 2 n 2 



2g s <f)(t E ) + MV t B 



1/2 



For M 2 — >• we define 



with 



1(0,8) =:£ 



/(") 



n=l 



2.s 



(56) 



(57) 



and 



/(n) = —27m 2 Inn + #(n) 



2 /4tt 2 2g^tE)\ 



+47rn 



47r 2 n 2 



+ 2^ s 0(ti 



1/2 



J E ■ 



(5? 



(59) 



The function J(0, s) can be evaluated using the Abel-Plana formula [331 132] 



7(0,5) 



E 27m 2 Inn Z" 1 , 
— - — — + / dng{n) 



n=l 



n 



dn 



9(.n) 



n 



2 s 



dy 

'o 



■,2-Ky 



-g(0) + O(s). 



(60) 



The sum on the right-hand side of (160]) will not affect the ^-function, since 
it corresponds to an M-independent term of 7(M 2 ,s). We have retained in 
( 160]) the regularizing factor 1 / n 2s only in the sum and in the second integral, 
since the remaining terms are finite for s — > 0. The splitting of the integral at 
x = 1 is made for convenience and does not affect the final result. In order 
to regularize the remaining integral, we have to integrate by parts several 
times. For this purpose we change the integration variable to x — 1/n and 
define the function 



F(x) = x 3 g(l/x) 



(61) 
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which is analytic at x — 0. The function 1(0, s) can be expanded according 
to 



7(0, s) 



I H (0) + 0(s) . 



(62) 



Performing integrations by parts and discarding the boundary terms where 
we cannot interchange the limits x — > and s — > [32], we find 



and 



7 R (0) 



J polc (0) 



25 d 4 F{x) 



1 d 4 F(x) 



288 dx 4 



dy 



x=0 



48 dx A 



1 

24 Jo 



(63) 



x=0 



In a; 



d 5 F(x) 



oo 

5«o) - £ 



27m 2 Inn 



n=l 



n 



2.s 



(64) 



The expression (163 p and the second integral in (I64p result from the convergent 
integral that remains after the integrations by parts. A factor x 2s /s in this 
integral can be expanded according to x 2s / s = l/s + 2\nx + 0(s) which leads 
to the aforementioned terms. 

The first term in fl64D results from the finite contributions of the boundary 
terms where the limit x — > and s — > can be interchanged. We find the 
explicit expressions 



J pole (0) 



g 2 s LH E p(t E ) 



and 
I R (0) 



24tt 2 

+ g B (f>(t E )t E ( 14 + 3 In 



- 23? \ iix 



+ 8tt 



2y/2ng a (l>(t E ) 

L 2 g s <p{t E ) 



(65) 



(66) 



2% 2 



2hg 2 s L^ 2 {t E )t E 
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Furthermore, we need the pole part of I(M 2 ,0) for M/0. Equation (ISTij) 
also holds for arbitrary M if we replace 4ir 2 /L 2 — > 4ir 2 /L 2 + M 2 in the 
definition of g(n). More concretely, we may expand I(M 2 , s) for large n and 
use the fact that the Riemann ^-function has a pole at s — 1, 

CR(2a + 1) = 1 + O(a ) ■ (67) 

We find 

JP ole (M2H _ 

(f^ + M 2 ) 7 

In the limit M — >• oo we find the regular part of i#(M 2 , 0) to be 

oo) = -f) (69) 

n=l 

This sum appears also in 1(0, s) and will therefore not affect £(0) and C'(0)- 
The logarithmic contribution of J(M 2 ,0) vanishes, since Cr(— 2) = 0. 

Using (H9|) and Eq. fll31j) from the Appendix, we can now compute the 
modified tunneling amplitude from J(M 2 , s). Expanding (I66p up to the first 
order in g s and using the normalization defined through (j55p . we find the 
modified tunneling amplitude 

T = r exp L^(tE)L 2 + gM^tA {7Q) 

Let us illustrate this result by using a definite expression for the function 
4>{t E ). With the choice 

W B ) = ^± tanh(c^) + *=±*±, (71) 

where 1/u is the characteristic time of the instanton, fTTO]) becomes in the 
limit of large Euclidean nucleation time T , that is, T ^> 1/cj, 

r-r exp^-^— + . (72) 

The first term in the exponent, which results from subexponential terms of 
the WKB-expansion, is negligible in the limit T 3> L. 

The result (j72|) deserves some explanation. The appearance of the quan- 
tization length L is due to the fact that the environmental field enters the 
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interaction quadratically. Integrating over all momenta leads to an expression 
which increases with the quantization volume. If we had chosen a bilinear 
coupling, the end result would have been independent of the quantization 
volume, since the environmental field would fluctuate around zero. 

If we neglect the first term in the exponent of (172|) . the tunneling ampli- 
tude will always be enhanced, since the product g s (p is positive; a negative 
sign would render the model unstable. 

Recalling the results of [15] , one would expect a reduction of the tunneling 
rate due to dissipative effects. In [15] the interaction with the environment 
leads to an effective friction term in the equation of motion for the system 
variable. Since we consider an interaction that is quadratic in the environ- 
mental field, the situation here, however, is differently] 

A finite number of environmental degrees of freedom (N harmonic oscilla- 
tors) would lead to a suppression of the tunneling rate, that is, to a negative 
sign in the second term of the exponential in (172|) . The mathematical reason 
is that for N oscillators we have in the leading term the sum 

N N 

5> = -(;v + i), ( 73 ) 

71=1 

which is always positive. In contrast, an infinite number of environmental 
degrees of freedom leads to a sum that must be regularized, e.g. with the 
help of the Riemann ^-function as above. One then gets 

n=l 

that leads to the positive sign in ( T72]) . Physically this can be interpreted 
analogously to the Casimir effect: due to the boundary conditions defining 
the functional determinant, there are fewer environmental modes than there 
would be without the boundary conditions, and thus there is less decoherence. 
This is a result that one would not expect on purely quantum mechanical 
(as opposed to field theoretical) calculations. 

So far we have restricted ourselves to tunneling in the flat Minkowski 
background. In a FRW universe with scale factor a, the eigenvalue equation 
(15TS|) changes to 

( d 2 3d d 4n 2 n 2 „ , . A , , . , , , N 

--772- - — IT + -TTY + 2 9*4>{t E ) u -A, t E ) = A u -A, t E ) . 75 
\ dt E a dt E L z a z J 

In this connection, note also the recent paper |34) where an enhancement of the tun- 
neling rate was found for the same mathematical reason, though in a different physical 
situation, with the electromagnetic field instead of an environmental one. 
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With the ansatz u(—X,t E ) = (pi{—X,t E )if2{t E ) and an appropriate choice of 
4>2{tE), we eliminate the first derivative in (I75|) . 

d 2 4vr 2 n 2 9a 2 , , 3 d f a\\ . , 



dt% L 2 a 2 4a 2 ys ^ v *" 2dt E \a, , 
= X Vl (-X,t E ). (76) 

The scale factor is in general a complex function of the Euclidean time t E 
that leads to complex eigenvalues of the differential operator. Assuming that 
no eigenvalues lie on the negative real axis, we can still apply the regular- 
ization method. The functional determinant related to (j76p is obtained by 
performing the substitutions 

L 2 -»■ L 2 a 2 (77) 

and 

2 9>#E ) ^(t E ) + - + _(- j. (78) 

We restrict ourselves to flat de Sitter space with a/a = —iH = const and 
with a small Hubble parameter, that is, Tq,L <C 1/H and H uj. Using 
fl56|) . we find a correction term of order g s H 2 , 



g s L(fi_T g s L 2 



T exp 



V 72 128tt 32tt 2 V V £ 



12 


8tt 




3L 4 


72 


1287T 



(79) 



where we have omitted terms that can be neglected for large To. In order 
to discuss the result quantitatively, we depict in Fig. [1] the dependence of 
ln(r/ro) on the length L which is roughly the size of the nucleating vacuum 
bubble. The Hubble parameter is set to zero, and we have evaluated ln(r/r ) 
using the exact expressions (I6"6"|) and (|68[) . Obviously, the correction term of 
the nucleation rate increases with L and the coupling g s . In Fig. [2] we depict 
the correction term of the exponent in (179"]) for small values of the Hubble 
parameter, that is, L,Tq 1/H. The finite Hubble horizon leads to a 
reduction of the exponent for small L and T . We note that for unrealistically 
strong coupling one gets a suppression of the tunneling rate (not shown in 
the figures); this is reminiscent of the quantum Zeno effect P^JfTB]. However, 
in that limit is may no longer be realistic to neglect the back reaction on the 
background. 
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Figure 1: We set L = T for simplicity and show the change of the modified 
tunneling rate with increasing L for different couplings g s . (For this choice, 
the exponents in (172]) and (f79|) still remain positive.) 



4 Cosmic Landscape 

The cosmic landscape motivated by string theory was discussed in various 
publications [TOl [TTj . Usually one considers Coleman-De Luccia tunneling 
[3TL l35| l36| 1371 138] between a huge amount of vacua and discusses various 
solutions of ad hoc rate equations. Under certain circumstances a continuum 
limit of these rate equations can be derived [53, BP] . 

Furthermore, finite temperature effects have been considered in [2TJ based 
on Hawking-Moss tunneling [22] • Rapid tunneling was proposed under the 
assumption that resonance tunneling is dominant in the landscape [221 1221 
I2S], see also [26] for a critics based on standard quantum field theory. In the 
following, we want to extend our model discussed in the preceding sections 
to multi-level systems in order to see under which circumstances an ad hoc 
rate equation can be formulated. 

One can, for example, generalize the Hamiltonian (1 12]) in the following 
way: 

H<f, = #dia g + H A = diag(cui, ...,u n ) + ^ Ajj\i){j\ , (80) 

Him = g s a 3 (t)diag(S u ,---,S nn ) J d 3 pa(p)a(-p) . (81) 

The interpretation is as follows: The numbers Ui denote the different lo- 
cal vacua of a cosmic landscape, and the Ajj are tunneling matrix elements 
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Figure 2: We set L = T for simplicity and show the change of the modified 
tunneling rate with increasing L for different Hubble parameters H. The 
coupling was chosen to be g s = 0.1 



which can be computed in WKB-approximation. The entries Sa in the in- 
teraction Hamiltonian distinguish the different vacua, which is an obvious 
generalization of measuring the left and the right well in the double-well 
system discussed above. 

Since the tunneling matrix elements are usually exponentially small, the 
short-time dynamics (short with respect to the tunneling times A^ 1 ) is de- 
termined by the decoherence rates. The off-diagonal elements of the density 
matrix now read, cf. (1261). 



Pij (t) = Pij (0) exp \- y ' K " 4 33 } Tr \± > ^ 1 j - i^A . (82) 

One can conclude from this expression that the suppression of interference 
terms depends crucially on the distance between different minima in the 
landscape. 

If one neglects possible degeneracies and assumes that the typical deco- 
herence rate is much larger than the tunneling rate, the system dynamics is 
determined by the equations [13] 



Pu(t) = - [ d8[H A (t),[H A (8),p(8)]]ii. (83) 

Applying the Markov approximation, one obtains 

pa(t) = \^2\A ik \ 2 (p kk - pa) , (84) 
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where A is chosen such that the coarse-graining in time is not too small and 
the approximation is valid. 

There exists another generalization of the model where the tunneling 
between different vacua may, in fact, be mediated by the environment. This 
is described by setting Ay = and Sij ^ 0, i ^ j in the above Hamiltonian. 
Indirect coupling between metastable states is a well-known phenomena in 
glasses, see, for example, [J7]. 

In the following we shall derive the master equations for the environment- 
mediated tunneling. The interaction Hamiltonian in the interaction picture 
has the form 



J5) 



Restricting ourselves to flat slices through de Sitter space, the operators 
cr(p, t) are given by 

a(p,t) = / p (t)a p e ipx + h.c. (86) 

with 



Mt) = ^7^W) +tH]e '^- (87) 



Applying the Redfield approximation [T3], we find for the system density 
matrix the expression 



^(t) = -Tr B / ds[HL(t),[HL(s),Ps(t)PB]]. 



In the limit of vanishing temperature, the bath density matrix is just ps = 
|0)(0|. The coefficients of the density matrix satisfy the system of differential 
equations 



w - ■ ■ — W - ■ 

ii j 3 1133 



39) 



Pji = 6 a Pkk( w iuk + w kuk) - Pji 

k^i 

with the correlation functions 

™timn(t) = 9s f dse^-^^S kl S mn a*{t)a\s) 
Jo 

x Jd 3 pJ d 3 q(a(p,s)a(-p,s)a(ci,t)a(-q,t)) (90) 
24 



and 



^nki(t) = 9 2 s f dse^^S H S m y(t)a*{s) 
Jo 

x J d 3 p J d 3 q(a(p,t)a(-p,t)a(ci,s)a(-q,s)) . (91) 

In deriving (189|) . several approximations have been performed: the Born 
approximation, which states that the total density matrix can be written 
approximately as a tensor product of the bath density matrix and the system 
density matrix, and the rotating wave approximation, which is valid when 
the intrinsic time scale of the system is much larger than the relaxation time 
of the open system. Since the correlation functions are not homogeneous in 
time due to the scale factor, the master equation is not Markovian. 

The rates in 089)) need not be exponentially small and may therefore 
dominate the dynamics of the string landscape. The transition probabilities 
between the vacua are symmetric, since we assume that the environment is 
described by a Gaussian wave function rather than an ensemble of states. 
On the other hand, the tunneling rates in the Pauli equations (1891) are not 
symmetric and jumping to lower energy levels is more probable than jumping 
to higher energy levels, depending on the bath temperature. It would be 
interesting to see how the situation changes if the Gaussian is replaced by a 
(micro) canonical ensemble. 

Evaluating the correlation functions we find 

< mn (t)=9 2 s f dse^-^ s ^S kl S mn a\t)a\s) 
Jo 

J d 3 k(^(f k (t)f:(s)r + j^y 3 J d* P \f k (t)\y p (s)\^ (92) 

and 

=fl? f dse^-^ s ^S kl S mn a'{t)a^s) 
Jo 

X ^V~ I ^( 2 ^( s )A ! W) 2 +(^p / d 3 P \f k (s)\ 2 \f p (t)\ 2 ) -(93) 

The dominating contributions in the correlators (1921) and fl93|) are given by 
infrared contributions k < Ha, since the phases in the integrands are oscil- 
lating rapidly if k > Ha. Neglecting the second term in the parenthesis of 
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f[9"2"j) . respectively (JHSJ, and applying the approximation e tk / Ha ~ 1 leads to 



t»LW = / dse^- w «><-*) 5 fc ^ mn a 3 (t)a 3 (s) (94) 



o 

X (^/ A (2^(4 + !ff ) 2 (^- !ff ) 2 <95) 

and an analogous expression for w~ nW (t). The evolution equation for the 
off-diagonal elements read in the Schrodinger picture 

/f h = -(Iy - i(uH - Wy ))pg* (96) 



with 

r y = 3? [ . + u£ K ) - «;+ i - ^ ] . (97) 



The imaginary part of the correlation function has been absorbed into the 
frequencies Wj. For = SaSy we obtain for a(t) ^> a(0) 

ry „ 9<2(5 ,,_^Z^)M _ M, (98) 



9(27r) 2 V^min ^m; 

where /c min is an infrared cutoff and fc max ~ if a. Solving fl82|) using the rates 
gives for large times 



\Pij(t)\ = Plj (0)exp[-g 2 s / dt'T i3 (t') 



to 



iimi 



This result can be compared with the off-diagonal element (1261) . Using the 
dominant contribution of the exponent given by ( 136ft . we find 

( ?(<!>+ - (j)-) 2 Va 6 H 2 ( 1 1 \\ , v 

[p ± (t)l=p ± (0)exp -g s ^ + ^ -3 — )) , (100) 

\ uu " \.rmm .rmax/ / 

which coincides with (199]) if one identifies (Sjj) with 0+ (0_). 

The transition between the different localized vacuum states is given by 
the rate equation 

pf? = -2 J2^<ktPT ~ ^tuuPlt) • (101) 



VI 
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For a(t) > a(0) we find 



X™ kiik - - (27r)2 9R2 + {ui _ Uk)2 k sJ ■ ^ 

Depending on the physical situation, that is, depending on whether a bath- 
induced coupling between different vacua or the tunneling dominates, (llOip 
or fl84l) describes the evolution of the cosmic landscape. 

This evolution is described by n coupled ordinary differential equations 
and can obviously not be solved for n ~ 10 500 vacua; therefore, further as- 
sumptions and approximations are necessary in order to obtain some physical 
insight. 

The approximation of Markov equations by Fokker-Planck equations is, 
for example, described in [49J and can always be applied if there is some 
small expansion parameter, for example, the ratio of the jumps between 
different vacua and the size of the tunneling landscape. A large number 
of vacua motivates the transition from discrete values pu(t) to a function 
Px(t), where x is a continuous coordinate in a smooth cosmic landscape. If 
the landscape is one-dimensional, one might consider the following scenario: 
There are probabilities to go to the left and to the right, which are described 
by functions a{x') and j3(x') if the observer is located at a position x'. These 
functions are the continuum limit next-neighbor transition rates in (}84l . 



respectively (jlOip . 



A|A M+1 | 2 = A^/3(x), 
A|A M _!| 2 = a>i -¥ a(x) , 



2Mwf i+l;i+1>i = (3 i ^(3(x), 



(103) 



(104) 



The probability that there is a local minimum in the potential between x and 
x + dx is Q r f(x)dx, where Q denotes the size of the cosmic landscape. There- 
fore, we take into account the tunneling rates and the distance to "nearest 
neighbors" of local vacua. 

The transition from the sum in (llOip to a continuous description can be 
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performed as follows: 



k k 




dx(5(x — kx) + ... + 5(x — k n ))w(i, x)p(x) 
dx J dkf(x — k)w(i, x)p(x) 
dxj(x)w(i, x)p(x) . 



(105) 



Following the treatment of [H] and assuming detailed balance, 
(3(x')'j(x)p(x') stat = a(x')^(x')p(x) stat , 



(106) 



where p(x) stat is some stationary distribution, the Fokker-Planck equation of 
diffusion type holds: 



This equation describes diffusion in an inhomogeneous medium, since the 
rates do not prefer a special direction in the cosmic landscape, that is, the 
dynamics is a random walk in an inhomogeneous medium. The drift term 
in f)107p is due to inhomogeneities in the cosmic landscape and vanishes for 
7(x) = const. In the following we will rescale the time such that the factor 
2/fl is absorbed. If the rates are time-dependent as in (11011) . the diffusion 
equation acquires a time-dependent factor on the right-hand side. 

Let us illustrate the solution of fll07p in two simple examples. If the pure 
tunneling given by (1841) dominates and 7(x) = AA 2 , where A is a typical 
tunneling rate, we obtain the usual solution for the diffusion equation, 



If the dynamics is environment-induced and given by f llOip . the diffusion 
depends on the scale factor. For a scalar-field environment and assuming 
a(t) = exp(Ht), the result is approximately given by 



dp(x,t) _ 2 d 1 d p(x,t) 
dt Qdxj(x) dx p s tat(z) ' 



(107) 




(108) 




(109) 



with 



D 



H 2 



(110) 
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where S denotes a typical transition element in the interaction (|8ip . 
Therefore the diffusion process may become faster due to the growth of the 
scale factor. This is, of course, only possible if the cosmic landscape and its 
environment exchange enough energy to lift the scalar field from one local 
minima to another. 

We have shown in our paper how a small positive value of the cosmo- 
logical constant could be justified through decoherence. A realistic scenario 
enabling the calculation of the observed value can, however, only be pre- 
sented after a definite cosmic landscape for the potential has been retrieved 
from a fundamental theory. 

5 Conclusion and Outlook 

In our paper we have addressed the problem why the cosmological con- 
stant (dark energy) has a small positive value instead of being exactly zero. 
Yokoyama had suggested that this could be due to the dark-energy field not 
being in its ground state (whose energy is assumed to be zero) but in a 
localized state [9]. Using a quantum mechanical model with a double- well 
potential, we have justified the localization of the dark-energy field in one of 
the minima. The crucial mechanism for this is decoherence - the emergence 
of classical properties (here, a localized state) by irreversible and ubiqui- 
tous interaction with irrelevant degrees of freedom ("environment") leading 
to quantum entanglement between the dark-energy field and environment 
[T2] . This localization is similar to the emergence of a chiral state for sugar 
molecules. 

More precisely, we have considered a Yukawa interaction between the 
dark-energy scalar field (quintessence) and environmental modes that can 
be interpreted as gravitational degrees of freedom. This interaction induces 
a dynamical suppression of interference terms connecting the two minima. 
Consequently, the ground state (which is a superposition of the two localized 
states) will for the dark-energy field evolve into an ensemble of localized 
states with an effective energy greater than zero. The decoherence factor 
is dominated by low-frequency super-Hubble modes. These modes become 
strongly entangled with the dark-energy field because their effective coupling 
involves an additional factor of Ha/p compared to sub-Hubble modes. 

Motivated by recent ideas in string theory, we have generalized our model 
and have considered an arbitrary number of perturbative vacuum states. 
Again, strong decoherence leads to the localization in a particular potential 
well. But one can also obtain the interesting limit of an environment- induced 
transition between different minima. Within the Markov approximation, we 
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have derived rate equations describing the latter feature. These transitions 
can dominate the time evolution if the coupling is strong enough to transfer 
sufficient energy from the environment to the scalar field in order to lift the 
latter out of the local minima. In the continuum limit we have found the 
Fokker-Planck equation for the distribution of the dark energy within the 
landscape of possible vacuum states. 

Another issue addressed in our paper is the change of the tunneling rate 
induced by the coupling to the environment. Is is known from various models 
involving bilinear couplings between system and environment, for example 
from the Caldeira-Leggett model [15] , that decoherence can lead to the sta- 
bilization of metastable states. This is analogous to the quantum Zeno effect, 
that is, the prevention of a decay process due to the continuous monitoring 
by the environment. 

We have, however, considered in our model a tri- instead of a bilinear 
coupling and have used a field-theoretical renormalization procedure instead 
of an ad hoc chosen cutoff frequency. For vacuum bubbles smaller or roughly 
equal to the Euclidean nucleation time, that is, L < 2irT Q /3, we have found an 
enhancement of the nucleation rate instead of the usual suppression, whereas 
for L > 27rTo/3 we have found a suppression. The enhancement is a conse- 
quence of field theory and similar to the Casimir effect. 

Quantum interaction with the environment gives a natural reason for the 
localization of the dark-energy field in a potential well and, consequently, for 
a small positive value of the cosmological constant. As long as the precise 
form of the potential is not known, its exact value can, however, not be 
computed. 

The present work can be extended in various directions. For example, it 
would be interesting to study the localization process when the dark-energy 
field is spatially inhomogeneous, that is, when it adopts different values in 
distinct spatial regions. An additional difficulty would then be the inclusion 
of collisions between different vacuum bubbles. 

Another possible extension would be the application of these ideas to the 
inflationary era in the early Universe, where the energy scale is much higher. 
The principal mechanisms of localization and modification of tunneling rates 
are the same, but the quantitative details are different. From these details one 
should be able to learn whether or how decoherence can lead to a metastable 
de Sitter solution and how this could affect the dynamics of the transition to 
non-inflationary eras following inflation including (re) heating. We hope to 
return to some of these issues in future publications. 
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Appendix 

In this section we shall give a short review of the ^-function renormalization 
method as presented in [32]. The authors there considered the renormal- 
ization of a functional determinant defined by a second-order differential 
equation. Usually, neither the eigenvalues nor the eigenf unctions of the dif- 
ferential operator are known exactly. Moreover, even if all the eigenvalues 
are known, the determinant is an infinite product of eigenvalues, which is in 
general a divergent quantity. 

To solve these problems, one represents the functional determinant via a 
generalized Riemann (^-function. The determinant of an arbitrary differential 
operator T> can be written as 

(Det£>) 1/2 =exp ^lnfJ A j = ex P^X^ lnA J = exp(W) , (111) 

where the eigenvalues of the operator are denoted by A. We define the gen- 
eralized ^-function through 

A 

which is a convergent series for some s > and can be continued analytically 
to s = 0. The exponent W in equation (11 IIP can be obtained through 

(113) 

s=0 

Since the eigenvalues of T> have the dimension of mass squared (recall that 
h — 1), this leads to a wrong dimensionality for W . Therefore we have to 
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replace ( I113p by 



1 d ^ f^ s 



A 



= -^C / (0)-^C(0)ln/i 2 , (114) 

s=0 



where we have introduced a renormalization parameter /i with mass dimen- 
sion one. 

The differential operator corresponding to a single field mode may be 
labelled by n. In quantum mechanics we are confronted with a finite number 
of modes, whereas in field theory we have to deal with an infinite number. 
For each fixed n, the eigenvalue equation reads 

D n u n (-X,t) = Xu n (-X,t) , (115) 

where A is determined by the boundary condition 

u n (-A,t o )=0. (116) 

This boundary condition together with a normalization determines the eigen- 
functions uniquely. All the boundary conditions (1116P can be collected in the 
equation 

Detu n (-A,t ) =0 , (117) 

where the determinant is taken with respect to all modes n and all eigenvalues 
A. With the help of the Cauchy formula, the generalized ^-function can be 
expressed as 

^=^IMV Du ^ toh (n8) 

with the contour C encircling all roots of equation (11171) . Deforming the 
contour C to a contour C which encircles the branch cut of the function z~ s , 
we find 

sin(7rs) f°° dM 2 d ^ , . - , 

We first consider the regularization method for a quantum mechanical sys- 
tem. The necessary information for the regularization of a system with a 
finite number of modes is contained in the function 

J(M 2 ) =J2 In MM 2 , t ). (120) 
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Expanding this function for large M leads to 

N 



I{M 2 -> oo) = ^(4 + I k In M 2 )M 2fc + (J R ) ln In M 2 + / R (oo) , (121) 



k=l 



where (I R )i n is the coefficient of the logarithmic asymptotic term of J(M 2 ), 
and I R (oo) is the asymptotic value of the regular part of J(M 2 ). According 
to [22], the (^-function can be expanded as 

C(s) = (I% n + s [I^ + 0(s 2 ), (122) 

where [/ R ]§° = 7 R (oo) — J R (0). As a demonstration we will apply this method 
to the harmonic oscillator. The eigenvalue equation 

d 2 \ 

■— + u 2 \ u(-X,t) = Xu(-X,t), u(-X,t o ) = 0, (123) 

has a solution of the form 

u(-X,t) = Ae V ^ Xt + Be~ V ^ t . (124) 

Performing the analytical continuation to the complex plane, A — > z, the 
function u adopts on the negative real axis the form 



u{M 2 ,t) = Ae Vuj2+M2t , (125) 

where we have neglected the exponentially decreasing term. Using for con- 
venience the normalization it'(O) = 1 leads to 



I(M 2 ) = -- \n(oo 2 + M 2 ) + Vu 2 + MH . (126) 



Since the term proportional to exp(— \/u 2 + MH) has been neglected, the 
analytically continued eigenfunctions do not respect the boundary condition 
u(-A,t o ) = 0. 

From (I126P we find J R (0) = — In a; + cut, (/ R )i n = —1/2. Using equations 
{HQ), (HH, and (1122)1 we arrive at 



Det 



d 2 2 
■df 2+U 



-1/2 

^e-f (127) 



which gives for large t the correct ground-state energy cu/2 for a harmonic 
oscillator, cf. Eq. (2.16) in the second reference of [3T] . 
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In general, the exact shape of the eigenfunctions u n (—X,t) is unknown 
and one approximates the u n with a uniform asymptotic WKB-expansion. 
This asymptotic expansion has the property that 

In u n (M\ t Q ) = WKB (n 2 , M 2 /n 2 ) (128) 

is uniform for M 2 /n 2 — > oo and M 2 /n 2 — > 0. In addition, it is also possible 
to use (I128p for the regularization of functional determinants in field theory, 
that is, if the mode number n is not bounded. The expansion (I128p has at 
most a finite power-law order growth in n [321 ED]- This fact allows us to use 
the parameter s to cure the divergences arising from the infinite number of 
modes. Changing the integration variable from M 2 — > n 2 M 2 leads to 

with 

/(M 2 , S ) = ^-^ln Mn (MV,t ). (130) 

n 

For a finite parameter s > the expression (11291) is finite. Analytic con- 
tinuation of the ^-function from its convergence domain to s = leads to 



COO = ^(/ P ° 1C )ln+(/ R )ln+[/ POl T 

+s! [ [ix - f dM2 ^ M2 dIP °dM^ 2) } + ° {s2) ■ (131) 

The coefficients (J polc )i n , J pole (oo) and J R (oo) are defined through the large 
M-expansion 



s 

+ J R (oo) + (J R ) ln In M 2 + 0{M 2 ) , (132) 
and the pole part 7 pole (M 2 ) is defined through 

I(M 2 ,s) = i L + 0(s ). (133) 

s 

J polc (0) and J R (0) are determined by 

/(M 2 ->0,s) = ^ + J R (0) + C(s). (134) 
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It is also possible to apply the regularization method if the differential 
equation exhibits singular coefficients. According to Olver |SD], the WKB 
expansion of a second-order differential equation 

d 2 v(M 2 t) 

{ df2 1 ] = uj{tfu{M\ t) = [f(t) + g(t)]u(M 2 , t) (135) 

has the form 

ft 



u(M 2 ,t) = C(M)(^(t))- 1 / 4 exp | J dt'{g{t')fl 2 j [I + h{t)\ (136) 
with h(t) = 0(M~ l ). The function h(t) can be expressed as Volterra integral 
h(t) = ^Yl-ex P |2^%V2 (n ^_2^^ ( ^^|j x 

xi>(t')[l + h(t')]dt' (137) 

with 

m ~9^(t) gV*(t) dt 2 g^{t) ■ 

Therefore the WKB expansion only makes sense if the kernel of (11371) is 
bounded. This leads to the condition 

q?(t) = [ dt'\ip(t')\ < oo . (139) 

Here, to an d t are the boundaries of the interval under consideration. The 
split of uj(t) 2 , see (If 35|) . is chosen such that singular coefficients like 1/t in 
the differential equation do not destroy the WKB expansion. This is the 
reason for the f/4-trick used in [32] . 

As already mentioned above, the evaluation of I(M 2 ,s) involves an im- 
portant approximation. In order to fulfill the boundary condition (If f 61) . 
two linearly independent solutions of the corresponding differential equation 
are required. After analytical continuation, the solutions are of the form 
u ~ exp(Mt) and u ~ exp(— Mt). The second solution is exponentially de- 
creasing and will therefore be discarded. This implies that the analytically 
continued functions u(M 2 ,t ) do not respect the boundary condition (II 16)) . 
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